Measurement of two-mode squeezing with photon number resolving 

multi-pixel detectors 
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The measurement of the two-mode squeezed vacuum generated in an optical parametric amplifier (OPA) 
was performed with photon number resolving Multi-Pixel Photon Counters (MPPCs). Implementation of the 
MPPCs allows for the observation of noise reduction in a broad dynamic range of the OPA gain, which is 
inaccessible with standard single photon avalanche photodctcctors. © 2013 Optical Society of America 
OCIS codes: 030.5260,040.5160,270.6570 
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Multiphoton entangled states attract considerable at- 
tention in modern quantum optics, as they represent 
promising resources for quantum communications, quan- 
tum computing, quantum lithography, and quantum 
imaging [1,2]. One of the most accessible ways to gener- 
ate such states is by using parametric down conversion 
(PDC), a non-linear process producing pairs of strongly 
correlated photons. The state produced by frequency 
non-degenerate PDC is a two-mode squeezed vacuum 
(SV) state written as follows [3]: 
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where \n)j denotes the Fock-state of n photons in the 
j-th mode, j = s,i denote signal and idler modes re- 
spectively, and Cn is the probability amplitude. Strong 
correlation between the photon numbers in the signal 
and idler modes, referred to as two-mode squeezing, re- 
sults in the suppression of the variance of their difference 
below the shot-noise limit [4,5] and can be quantified by 
the noise reduction factor (NRF) given by 

Var(iV, - N,) 
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where Ns and Ni are the photocount numbers of the de- 
tectors in the signal and idler modes respectively, and 
(O) = tr(/50) defines the mean value of an observable O 
for a given state p. The two-mode squeezing can be stud- 
ied with various types of photodetectors, depending on 
the photon fluxes of the PDC. At a moderate gain, when 
the SV state contains several photons per pulse, it is es- 
sential that the detector is able to resolve several simul- 
taneously impinging photons. Such detectors are referred 
to as photon number resolving detectors (PNRDs). 

Development of the PNRDs attracts a considerable 
interest in modern optical engineering. Several PNRD 
technologies are available up to date, including cryo- 
genic detectors [6-8], branched single photon avalanche 



detectors (SPADs) [9,10], hybrid photodetectors [11], in- 
tensified and electron- multiplied CCDs [12]. 

Here we consider a multi-pixel photon counter 
(MPPC) which is a commercially accessible PNRD [13]. 
In a MPPC, several hundreds of SPADs (pixels) are em- 
bedded into a chip which is illuminated by a diffused 
light beam, such that the chance of more than one pho- 
ton hitting the same pixel is negligible. Upon photon 
detection, a pixel produces an electrical signal, which is 
summed with the signals from other pixels at an out- 
put circuit, thus enabling the photon number resolution. 
Despite its attractive features, such as ease of opera- 
tion, low cost, and high quantum efficiency (q.e.), the 
crucial drawback of the MPPC is the optical crosstalk 
between pixels. During an avalanche in a pixel, some 
spurious photons are emitted, and eventually registered 
by the neighboring pixels [14]. As the crosstalk happens 
almost simultaneously with the detection of real pho- 
tons, its contribution cannot be distinguished from the 
actual counts, and therefore requires an accurate theo- 
retical modeling [14-16]. 

The performance of the MPPC has been studied in 
earlier experiments with coherent and pseudo-thermal 
light [14-18]. In [19] the MPPC was used in the measure- 
ments of the intensity correlation function of the SV 
state, however restricted to the regime of much less than 
one photon per pulse. In this work we expand the appli- 
cation of the MPPC to the study of a relatively bright 
SV state with up to almost 5 photons per pulse by the 
direct observation of the two-mode squeezing. 

For modeling a realistic MPPC, we follow [15], where 
loss and crosstalk are considered as Bayesian processes. 
Let us denote (n) as the mean photon number imping- 
ing on the detector, P as the crosstalk probability, and 
77 as the q.e. of a single SPAD pixel, which also accounts 
for the losses in the optical system. Here, beams of rel- 
atively high intensities have been considered hence the 
saturation of the MPPCs has to be taken into account by 
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assuming that each MPPC can only resolve up to A'max 
photocounts. We assume that the two MPPCs used in 
joint detection have the same values of rj, P, and iVmax- 
Then the positive-operator valued measure for detecting 
Nj photocounts in the j-th mode is given by 
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The variance of the difference of photocount numbers is 
given by [20]; 



Var(iV, - N,) 



m - {N^f + {ND - (5) 



where the p-moment operator of the photocount number 
operator is given by iVj = Y.'n'^o ^j^N, ■ Using (2) - 
(5), the NRF can be calculated for any given two- mode 
state. As (n) the model yields NRF = for 
the coherent state which is greater than unity-the value 
expected for a lossy PNRD without the crosstalk. Thus, 
the effect of the crosstalk is to increase the NRF. Also, 
as (n) — for the SV state the model yields NRF = 
- (1 + P)'n- Thus the difference in NRF values for 
the two states is equal to the effective q.e. of the detectors 
r]E = {I + P)r], differing from the absolute q.e. due to 
the presence of crosstalk [21,22]. Numerical calculations 
show that the effect of saturation is to decrease the NRF 
with increasing (n) (see Fig. 2). 

In the experiment the 4-th harmonic of Nd:YAG 
pulsed laser at 266 nm (Crystalaser, repetition rate 20 
kHz, 30 ns pulse width) was used as a pump for a travel- 
ing wave OPO (Fig.l). It was sent through two 5 mm 
length BBO crystals with oppositely directed optical 
axes in order to compensate for the spatial walk-off of the 
pump. The BBOs were cut for coUinear frequency non- 
degenerate type-I PDC with signal and idler wavelengths 
As = 500 nm and = 568 nm, respectively. A half-wave 
plate and a polarization beamsplitter (not shown) were 
used to control the pump power. After the BBO crystals, 
a UV-mirror filtered out the pump beam, while the PDC 
beam was transmitted. The signal and idler beams were 
split by a dichroic mirror; the former was reficctcd and 
the latter was transmitted. In each arm there were lenses 
with their front foci at the BBOs. Two iris apertures were 
set according to ds/di = A^/Ai in the back focal planes 
of the lenses, selecting the conjugated spatial modes [5]. 
The spectral bandwidth of PDC was calculated to be 3 
nm, and was restricted by the spatial selection. The PDC 
radiation was then collected by lenses and formed a beam 



spot of 1.2 mm (at level l/e^) at the MPPC chip. The 
ambient light was suppressed by means of interference 
filters centered at 500 nm and 568 nm (15 nm FWHM, 
peak transmission 95%, not shown). A 0.1 dB neutral 
density filter was introduced in the signal arm to bal- 
ance the q.e. of the MPPCs to the value of the q.e. at 568 
nm. Two MPPC modules (Hamamatsu, model 10751-02, 
20x20 pixels each of active area 50 x 50 /irn^, total chip 
size 1.5 X 1.5 mm?), were sealed into a custom housing 
with an AR-coated optical window and were thermoelec- 
trically cooled to -4.5°C to reduce the dark noise. The 
signals from the MPPCs were digitized by an AD card 
(National Instruments, PCT5154) and discriminated by 
their amplitudes. The AD card was synchronized with 
the pump laser with the detection window set to 70 ns. 
The data analysis was done with Lab View, and the theo- 
retical fits with Mathematica. The MPPC dark noise and 
ambient light of the setup were measured with an extin- 
guished PDC, for which the polarization of the pump 
was 90 degrees rotated by the half-wave plate, and then 
subtracted from the data. 
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Fig. 1. (Color online) The experimental setup: Nd:YAG laser 
at 266 nm pumps two BBO crystals. OA, optical axis; UVM, 
UV-mirror; DM, dichroic mirror; M, mirror; LI (f=500 mm), 
L2 (f=200 mm), lenses; Dl {ds = 10 mm), D2 (d» = 11.4 mm), 
iris apertures. The inset shows the setup for measurements 
of the coherent state: cw Nd:YAG laser at 532 nm is chopped 
by an acousto-optical modulator (AOM); ND, neutral density 
filter; NPBS, non-polarizing beamsplitter. 

Similar measurements were performed for the coherent 
light obtained by attenuating the 2-nd harmonic of a 
continuous wave (cw) NdiYAG laser at 532 nm (Photop, 
power 50 mW) and chopping it by an acousto-optical 
modulator at a frequency of 20 kHz with a pulse duration 
of 30 ns. 

Fig. 2 shows the experimental data and theoretical 
fits of NRF versus {fi) for the coherent and SV states. 
From the difference of the measured NRFs for the two 
states close to (n) = one finds the experimental value 
of the effective q.e. for the squeezed state at 568 nm. 
The value yields ?7^gv = O-l''' i 0.03, and it was fur- 
ther used for the conversion of photocounts to photon 
numbers. Effective q.e. for the coherent state at 532 nm 
was found from ??^gvj assuming the dependence of q.e. 
on the wavelength, provided by the manufacturer, and 
yielding t/^'^p^j^ = 0.20 ± 0.03. The theoretical fits of the 
data to the model yield iVmax = 3, P = 0.280±0.005, and 
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?7 = 0.163 ± 0.007 for the coherent state and TVmax = 3, 
P = 0.30 ± 0.02, and 77 = 0.145 ± 0.009 for the SV 
state. These parameters yield = '^■^l i 0.01 and 

r]^sv — 0-19 ± 0.01, which are in good agreement with 
the experimental values oi tje used for the calculation of 
photon numbers. The results clearly demonstrate that 
NRF gradually decreases with the increasing (n) for 
both studied states, which is caused by the saturation 
of the MPPC. At the same time, the NRF for the SV 
state lies below the one for the coherent state - a sig- 
nature of squeezing. The crosstalk elevates the NRF at 
low (n) to values above unity as predicted by the the- 
ory. Surprisingly, the found value of A^max appears to be 
much smaller than the one expected from the MPPC 
with 400 pixels. We have found that the MPPC modules 
do not allow reliable discrimination of photon numbers 
beyond (h) = 5 due to the limitations of the built-in 
electronic circuitry, which was inaccessible for any ad- 
justment. This might also be the cause of the deviation 
of the measured NRF from the theoretical prediction for 
(n) > 5. Note that reliable photon number resolution up 
to 20 photons was demonstrated with similar detectors 
using customized electronics [23]. 
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Fig. 2. (Color online) Dependence of the NRF on the mean 
number of photons measured for the coherent state (black 
circles) and the SV state (red squares). The theoretical fits 
for the NRF are done for the range of average photon numbers 
up to 5 photons (vertical dashed line), which marks the range 
of reliable photon-number resolution. The solid lines are fits 
to the theoretical model with = 0.9999 and = 0.9994 
for the coherent and SV states, respectively. 

In conclusion, we have performed an experiment in 
which the NRFs of coherent and SV states have been 
measured with commercially available MPPC modules 
at relatively high photon numbers (up to 5). The main 
conclusions we can make are: (1) the crosstalk in the 
MPPC leads to an increase in the NRF for both coher- 
ent and SV states, (2) saturation of the MPPC leads 
to the decrease of NRF with increasing photon num- 
bers. The experimental data agrees with the theoreti- 
cal model which takes into account the saturation and 



the crosstalk although the fits deviate at higher mean 
photon numbers, where the photon number resolution 
of the detectors is limited. These results extend the use 
of MPPCs for the characterization of quantum light in 
a significantly broader range than that of conventional 
SPADs. 

The work was supported by the A-STAR Investigator- 
ship grant. 
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